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The pion-pion scattering phase shift is computed using LapH propagators. The LapH method for 
computing quark propagators is used to form two-particle correlation functions with a number of 
different operators. Excited state energies of two-particle states on 2+1 dynamical, anisotropic 
lattices {m K ~ 390 MeV) are computed to determine the phase shift in the 1 = 2 channel. The 
signal for t-to-t diagrams for the 1 = channel are also presented to demonstrate the efficacy of 
the stochastic LapH method. 
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1. Introduction 

The determination of the spectrum of excited states in Lattice QCD with light dynamical 
quarks is hampered by the need to include explicit multiparticle operators in the variational ba- 
sis of interpolating operators. The need for finite momenta operators seems to require an all-to-all 
calculation of the quark propagator for the simulation of even the simplest two-pion state. However, 
it was pointed out [1] recently that hadron operators that are constructed from smeared propagators 
have a natural cutoff in momentum space and that this can be exploited to form one timeslice-to-all 
quark propagators if the quarks could be smeared. 

This method of computing the quark propagator below the cutoff defined by the quark field 
smearing ("distillation") was used to compute the Tin scattering length in the 1 = 2 channel on 
dynamical, anisotropic lattices and presented in Williamsburg. We extend this calculation by de- 
termining the excited states and the phase shift in this channel using the finite volume method in 
Euclidean space-time ([2]). Furthermore, we modify the distillation algorithm by combining it with 
a particular stochastic method known as noise partitioning ([3]) or noise dilution. This not only al- 
lows a more efficient way to compute all-timeslices-to-all propagators, but it reduces the linear 
volume dependence of the number of eigenmodes required to achieve the same quark smearing in 
larger volumes. This was a particular difficulty that had to be solved in the distillation method. 

We present here preliminary results for the 1 = 2 phase shift (see [1 1] for recent calculations) 
and the signal for the t-to-t pieces of the correlation function in the 1 = channel using this new, 
stochastic smearing algorithm. 

2. Construction of Operators/Correlators 

2.1 Quark Field Smearing 

We begin the construction of two-particle correlation functions with the low-mode filtering of 
the quark fields via the eigenvectors of the 3-dimensional Laplacian operator, A. We have computed 
the low eigenmodes (up to N ev = 128) of the three-dimensional, gauge covariant Laplacian operator 
on each timeslice, 

Av (i) = _A,vW 

where X$ is the eigenvalue with the smallest magnitude of the lattice Laplacian operator, 

M^) = t{ ^)S(y,x + k) + ul(x-k)8(y,x-k)-28(x,y)}. 
k=\ 1 J 

The tilde on the gauge fields indicate that they have been smeared using the stout-smearing algo- 
rithm ([4]). 

The quark fields are then smeared using the smearing operator, 

S to (x, x ') = ^&(-X i + a 2 )v^(to,x)®v^(to,x) 

i 

and the smeared quark field is given by \jf(to,x) = S to (x,x')y(to,x'). 
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We use these quark sources to compute the solution vectors of the Dirac matrix using standard 
methods. This requires solving for the solutions of each of the eigenvectors used to smear the 
source field, 

<j)®(t,x)=M- 1 (t,x;ta,x )v®(to,xo). 
The solution vectors, <j)( l '(t,x), are smeared using the same smearing method but on timeslice t, 

f i \t,x)=S t {x,x / )^ i \t,x / ). 

The smeared quark propagator is then given by, 

Q(t,x;t ,x Q ) = S t (x,x)M-\t,x;to,x )S to (x ,x Q ). 

3. 1=2 Phase Shift 

The quark propagator, Q(t,x;to,xo), is a ?o-to-all-? propagator. We can therefore construct pion 
operators with finite momenta without further inversions (i.e. at no extra cost). We compute the 
1 = 2 pion scattering phase shift by using several, finite momenta pion operators and diagonalize 
the correlation matrix to determine the energy eigenvalues of the two-particle state. 

3.1 Parameters 

We use the Nf = 2 + 1 anisotropic lattices with anisotropy tuned to a s /a t = 3.5 ([5], [6]). The 
lattice spacing in units of ro is given by ro/a s = 3.221(25) and the lattice size was 20 3 x 128. 
The results are from 90 configurations separated by 40 trajectories. A simple jackknife analysis 
suggests that the results are independent within the errors. Preliminary results for the two-pion 
correlation function and scattering length were reported on a smaller lattice, 16 3 x 128 with a pion 
mass of ~ 390 MeV in Ref. [7]. We use the same pion mass and five of the lowest momenta pion 
operators to project out the S-wave scattering state in the center of mass frame. 

3.2 1 = 2 Correlation Function 

The 1 = 2 correlation function is constructed in the usual way by computing the "direct" and 
"crossed" diagrams. The direct diagram is simply the square of the single pion correlation function, 

C n (t,to) = [Q(t,x;to,xo)] Q(t,x;t ,xo). 

The correlation matrix for the quark exchange diagram is given by 

cjf° M) (Mo)= (vl t M u l (z',t-y',t )V x , % ) {v;„ t V rtQ )e- k ^" (v_l t M u 1 (z' Az'" ,t Q )V z , % ) (V^) 

(vlM u \z,t-,x",t )V x% ) (y} t Vy t0 )e^y" (Vj ft? (x> ,t;/ ,t )V^ t0 ) (y} t Vy t ) 

where V x , t is the matrix whose columns are the eigenvectors, vW and the sum over the momenta 
is carried out to project out the S-wave. The 1 = 2 channel is obtained by subtracting the quark 
exchange diagram (d™ 1 )) from the square of the single pion correlation function, C n jC n j, 

Cij = Cx.iCjzJ — C}j (t , to) 
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Table 1: Phase shifts from the two-pion energies. 



n 


n 

u t r '3t3t 


a s p® 


h 


8(p®) 





0.150(3) 


0.04(4) 


0.015(33) 


-1.8(57) 


1 


0.275(4) 


0.326(4) 


1.078(26) 


-15(5) 


2 


0.368(4) 


0.457(6) 


2.118(58) 


-16(7) 


3 


0.447(9) 


0.54(1) 


3.01(13) 


-3(30) 



The correlation matrix is diagonalized at t* /a t = 25 with the metric timeslice att /a t = 15 so as 
to reduce the contribution from excited states as much as possible without losing the signal to the 
noise. The stability with respect to t* and t has been checked for the lowest-lying four states. 

3.3 Extracting the Phase Shift 

We follow Luscher's method to compute the phase shift in the infinite volume from the spectra 
of the two-particle state in a finite volume. First, the physical momenta of the pions is determined 
from the spectra of the two-particle states and the dispersion relation, 



(a s p„) 2 = £, 2 (a t m) 



( {a t En 



\ 2(a t m) 



1 



where E, = 3.5 is the anisotropy and a t m is the pion mass (in lattice units) at rest. We then compute 
the modified Zeta function, Zoo(l;n) to obtain the phase shift at momentum a s p, 



tan8(p„) 



71 



3/2 



where n = (a s p n ) 2 / 



Zoo(l;«) 

The results are tabulated in Table 1 and plotted in Fig. 1 



271 

L/a s 



4. t-to-t Diagrams 

One of the major challenges in lattice hadron spectroscopy is the evaluation of contributions 
from disconnected diagrams and box diagrams that appear in correlation functions. These diagrams 
require the quark propagator from a timeslice t, back to t {t-to-t diagrams) on some number of 
timeslices. One way to compute these contributions would be to compute wall propagators from 
every timeslice t which would require at least N t times more inversions than before (on our lattices, 
this factor is roughly 128). While this is not impossible to do in practice, we note that it may be 
unnecessary as the measurements from neighbouring timeslices may be strongly correlated, and 
also because the signal diminishes exponentially with time. The other more important reason to 
combine the LapH method with a stochastic algorithm is to control the number of eigenvectors of 
the Laplacian that needs to be computed as one progresses to larger and larger lattices. 

4.1 Stochastic Estimation 

We choose Z4 noise for our stochastic sources which fits naturally with having complex fields 
on our lattice. It has also been shown in some cases to have smaller variance than other noise 
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Figure 1: Preliminary results for the 1 — 2 phaseshift at m K = 390 MeV (Our points are shown as circles). 
The two-loop chiral perturbation theory curve was reproduced using the parametrization from Ref. [10] and 
the 390 curve was generated by shifting the pion mass. 



choices ([8]). Independent Z4 noise sources, p™, need to be chosen for each quark line in the 
correlation function of interest. These have the property 

«p»=0 and ((p [fl pj]» = ^ 

where the double bracket indicates an average over the noise sources. The quark propagator on a 
given configuration can then be written as, 

Q = DjS t M- l ((pp^))S t0 Dl 
= (DjS t M- l p (D k S !oP y)) 

where Dj is the covariant displacement operator in direction ij. 

In order to avoid contaminating the eigenvectors of the Laplacian with our random noise 
sources, we place the Z4 noise only in the LapH subspace. The noise vectors therefore only have 
spin and eigenvector indices on each timeslice. The noise vectors will be fully diluted in the time 
direction for connected diagrams, but can be interlaced in time for the disconnected diagrams. 
The full dilution scheme is given by 

P[A]«(0P[g]j'i'(^) = S A B8t t '8 ss '8w 

(without an average over the noise sources). However, this scheme is very expensive and is un- 
necessary in practice. The scheme dependence of some observables have been reported elsewhere 
in these proceedings ([9]). We only note here that the full time-dilution is usually unnecessary for 
computing disconnected diagrams since the noise on different timeslices do not interfere with each 
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Figure 2: Logarithmic ratios of correlation functions. The data with the label "1=0" has all contributions to 
the correlator except for the disconnected diagrams. 



other as long as they have a separation that is five (a conservative estimate) or greater. This leads to 
an interlaced dilution scheme which have noise sources on more than one timeslice, but separated 
by several timeslices. The same can be done for the eigenvector indices and the spin degrees of 
freedom. 

4.2 A box diagram example 

We have computed the box diagram that appear in the 1 = channel to test the efficacy of the 
"stochastic LapH" method. The box diagram can be written in a compact form by using the noise 
source vectors, Pms^oA), and the corresponding solutions, (f>\A] s ,i{t ,x) = M _1 Puh^o^o) as 

(/$!^))s*8S(*o) (/$Vi))s*g&)) (4]!>o)r4t^)) 

(4.1) 

where the colour indices are contracted within the round brackets and the dilution indices have 
been combined into one superscript. (This example is given for the usual Dirac matrix.) 

The logarithmic ratios for the 1 = correlation function with the box diagram but without 
the disconnected diagrams is shown in Fig. 2. We obtain a good signal for the box diagram with 
with a single timeslice source for the propagator from to to t (and interlace 6 for the eigenvectors) 
and time (interlace 12) and eigenvector (interlace 4) dilutions with full spin dilution. The ratio of 
correlators of neighbouring timeslices for the 1 = 2 and the 1 = correlator (without the completely 
disconnected diagrams) are shown are shown in Fig. 2. 

5. Summary 

The LapH quark smearing algorithm (a particular choice of distillation) has been tested for 
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two of the simplest two-particle states. For the 1 = 2 channel, the phase shift has been computed 
with several pion operators with non-zero momenta through Luscher's finite volume method. We 
have obtained a good signal for the first three, lowest lying momenta states in a (~ 2.5 fm) 3 volume 
for m n ~ 390 MeV with 90 configurations. It is clear that we require lighter pions and larger boxes 
in order to compare with chiral perturbation theory calculations. These simulations are under way. 

The box diagram in the 1 = channel has been computed using a new stochastic algorithm 
(stochastic LapH) in order to handle the t-to-t diagrams. Apart from the fact that this method can 
extend the LapH algorithm to large lattice volumes, it allows the simulation of t-to-t diagrams in an 
efficient way. We are currently working on the full 1 = calculation with disconnected diagrams 
and the 1=1 channel to study the p decay. 
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